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We present a general description of the propagation properties of quantum gravity modified elec- 
trodynamics characterized by constitutive relations up to second order in the correction parameter. 
The effective description corresponds to an electrodynamics in a dispersive and absorptive non-local 
medium, where the Green functions and the refraction indices can be explicitly calculated. The 
reality of the electromagnetic field together with the requirement of causal propagation in a given 
referrence frame leads to restrictions in the form of such refraction indices. In particular, absorption 
must be present in all cases and, contrary to the usual assumption, it is the dominant aspect in those 
effective models which exhibit linear effects in the correction parameter not related to birefringence. 
In such a situation absorption is linear while propagation is quadratical in the correction parameter. 
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CN ■ I. INTRODUCTION 

^ ' The different effective descriptions of quantum gravity effects at low energies can usually be expressed in terms of 
modified dispersion relations, with a polynomial dependence on the energy and the momentum. Such modifications 
00 include standard Lorentz invariance violations as well as possible extensions of Lorentz covarianceQ. In the case 
of electrodynamics these heuristic approaches can be interpreted in terms of an effective field theory describing the 
propagation of the electromagnetic field in an effective medium, provided that space-time coordinates and momenta 
^ ; commuted S S 1 S S • Otherwise, as it happens for example in double special relativity models, the ordering 
ambiguity in the Fourier transform imposes a different approach outside the scope of the present work|9j. In fact, most 
of the models considered in the literature, regardless of their quantum origin, are ultimately interpreted as a standard 
low energy effective theory when their proponents want to probe the corresponding phenomenological consequences 
by making use of the standard meaning of the required observables, such as energy and momentum. Of course, it is 
possible to abandon the effective field theory framework but then it is mandatory to revise the standard concept of 
■ observables together with unfolding the rules to measure them. 

In the effective field theory version of quantum gravity inspired modified electrodynamics the equations of motion 
O . acquire the structure of typical Maxwell equations in a medium lOj 

5^ ; ik B = Anp, k - B = 0, (1) 

W)' kxE-wB==0, ik X H + itjD = 47rj, (2) 

where the D and H fields are characterized by constitutive relations of the form 
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■ D' =a'^Ej+p'3Bj, (3) 

' W = P'^Bj+a'^Ej. (4) 

The coefficients a*-' , p*-' , /3*^ , and cr*^ are polynomials in lo and fc*. Once these constitutive relations are given, we 
can proceed as usual in the case of an electromagnetic field in a non-homogeneous medium. From a heuristic point of 
view, it is interesting to remark that these constitutive relations correspond to effective media that are non-local in 
space and time, which can be interpreted as a footprint of the granularity induced by quantum gravity. Although it 
seems that these media could be in principle quite arbitrary, the requirements of reality of the electromagnetic field 
together with causal propagation lead to constraints that must be satisfied by these effective theories on very general 
grounds. 

Our definition of causality at the effective theory level is the standard one, whereby the field strengths are obtained 
from the sources via retarded Green functions. It is important to remark that here we will deal only with the aspects 
of causality related to the dispersive character of quantum gravity effects as they are manifested in a given reference 
frame, a character related to a generalized susceptibility theorem. There is another more geometrical aspect related 
to Lorentz transformations connecting different reference frames, associated to the possible existence of closed causal 
curves which could appear when there are superluminal velocities. No matter how close to c the superluminal velocities 
are in a given reference frame, in a boosted enough frame these particles will propagate to the past. The analysis of 
the possible existence of closed causal curves will not be discussed here. 

In relation to Eqs. H1I2|) . we adhere to the point of view of Ref. (TJ where E, B, D, and H are considered as 
fundamentals fields, even in vacuum. In fact, the excitations D and H are the potentials for the charge density p 



2 



and the current density (j + dJD), respectively. Thus, they are directly related to charge conservation. On the other 
hand the electric and magnetic field strenghts E and B are the forces acting on unit charges. Our notation does not 
reflect their intrisic geometrical content described by the characterization of D and B as three-dimensional two-forms, 
together with E and H as one-forms. It refers to the corresponding components of the forms where the Hodge duality 
has been used to rewrite antisymmetric two-form components as one-form components. 

To introduce our effective field approach we will briefly recall the Lagrangian for the electromagnetic field in a 
non-dispersive medium. It is usually written as 

L = --^F.^x^^'^^^^f'^F^f, - j,A^, (5) 

lOTT 

where F^^i, = d^Ai, — d^Afj,, and x^t^'^^^"'^^ is a constant tensor that encodes the information about the medium. This 
structure, where the basic dynamical field is A^, warrants gauge invariance and hence charge conservation. The 
dynamics is given by the equations of motion 

a^ij^" = iTTf, (6) 

together with the constitutive relations 

^M-(^) ^^M[«/3]J7^^(2.). (7) 

Defining, as usual, the electric and magnetic field strengths as Fqi ~ Ei and Fij — --e.ijkBk respectively, and the 
corresponding components of the excitation H^" , H'^'^ — and H"^^ ~ —e^^^H^, the constitutive relations become 

= 2x["^l E, - x'"'! ^,nnjB, , (8) 

= e,;fcX^''""^'li?, - ^e»;fcx[''"™"lemn,B„ (9) 

and the effective equations of motion acquire the usual form for an electromagnetic field in a medium given in the 
relations H1I2(I . 

In the case of a non-local space-time medium this formalism is also suitable, but now we have 

L = --^[ d^ii^ (i-)^ (10) 

IbTT J 

instead of The constitutive relations become 

Demanding that the Lagrangian (|10|l be real for real fields F^fs (E, B) implies also the reality of H'"' (D, H). Writing 
-^[Mf][a/3] (^jfT — x"') in terms of its Fourier transform 

^M[«/31(2;'- - i'^) = y d^/fc g-^k-ix-x) ^M["/31(fc<-)^ (11) 

we can easily demonstrate that can also be written 

L = -^F.'^i^n (x^^-'i'^^k^^..) FMxn) ■ (12) 

This shows that for a non-local medium we can use a Lagrangian with the same form as (O, but now with 
being a derivative operator. Thus we are dealing with a higher order theory for the variables A^^. In particular this 
implies that dynamical functions, such as the density of energy and momentum, will depend not only on the fields E, 
B, D, and H, but also on their space-time derivatives. This will pose additional constraints regarding the causality 
of the system, since the causality of a given function does not warrant the causality of its space-time derivatives due 
to the presence of the energy cutoff in the effective theory (see Appendix A). 

If ig symmetric, in the sense that for each set of index values (n, a, (3) (no sum with respect to repeated 

indices) 

j d\ F^. (x'^'-l'^'^li^o/j) = I d\ F^p (x'"''l''^"li^p.) , (13) 
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is satisfied. Thus it is possible to perform integrations by parts rendering the equations of motion of the same form 
as the usual non-operator case. In the following we assume that this property holds. The above relation Ijl^l) can be 
rewritten as 

^ ;^Ml[M-](„,a^)^ (14) 

The expansion of the operator in terms of space-time components 

47r 47r Idtt 

can also be used to characterize a quantum gravity modified electromagnetic Lagrangian which preserves 3-d isotropy, 
provided that ^^I^^Ho*]^ T^[o»][ife]^ T^[ife][Oi]^ ^[ij][mn] ^^.g tensors under spatial rotations, but not the components of a 
y\tJ-v\[ap\ standard Lorentz tensor. 

In this paper we use this approach to discuss the main properties of quantum gravity induced effects in electro- 
dynamics, which are labelled by a parameter ^ usually considered to be proportional to the inverse Planck mass 
Mp. We provide a general description of such modified electrodynamics including expressions for the Green functions 
as well as the corresponding index of refractions, in terms of a parameterization of the constitutive relations up to 
second order in ^. Some preliminary work along these directions have already been presented in Ref. (l4l] . We also 
discuss the implications of reality and causality in the field propagation, applying them to some specific cases already 
considered in the literature. In a different context, causality in Lorentz invariance violating theories has previously 
been discussed in Ref. [l^ . 



II. CONSTITUTIVE RELATIONS AND EQUATIONS OF MOTION 

In the following we consider corrections up to second order in the parameter ^, assuming for simplicity that we are 
in a reference frame where there is invariance under rotations. This would correspond, for example, to the rest frame 
= (1, 0) in the Myers-Pospelov model 13]. Our approximation includes contributions from dimension five and six 
operators in the language of effective field theories. The modified Lagrangian yields the equations of motion: 



and the structure of the field tensor Fn,v leads to 



Comparing Eqs. H8I9|) and (|3I4|I we have 



5o 



f, - [iO] [Onx] ci I - [ iO] [mn] p 
9 V [*0] [O™] El _|_ ^ [ iO] [mn] p 



(16) 
(17) 

(18) 



X 



X 



^[Oi] [mn] 
\lk] [ran] _ 



(19) 
(20) 



In order to analyze the propagation of the fields and to define a refraction index it is important to make explicit the 
dependence of the constitutive relations on the frequency lo and the momentum k. To achieve this we expand their 
coefficients in space derivatives. Taking into account that these models can be understood as perturbative descriptions 
in terms of the parameter ^ = £./Mp we have, up to order 



(21) 
(22) 
(23) 
(24) 



where a^, I3ai cta, and pA, with A = 0,1, 2, are 5*0(3) scalar operators, according to the assumption that there is no 
preferred spatial direction. By relaxing this assumption it is possible to generalize the present approach to models 
that exhibit some preferred spatial directions. Examples of these last models are given in Refs. Ill, [11], and 
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where a Lorentz violating vector n'^ = (n'',n*) appears. Only in the reference frame where = do they yield an 
isotropic propagation, otherwise anisotropic effects occur. For simplicity, in the following we restrict our discussion 
to the isotropic case. 

The symmetry of x in Eq. implies that the terms in a*-' and /3'-' which contain an even number of derivatives 
are symmetric under i <-!■ j, while the terms with an odd number of derivatives are antisymmetric. In the case of p*-* 
and a'^ Eq. (O leads to 

= -CTj. (25) 

The property V-B = implies that P2 and p2 are irrelevant, and thus we take /32 = P2 = (and hence (T2 = 0). 
Thus Eqs. ii>lli>4|l finally reduce to 

a'^ = aov"^ + aiie'^^'dr + a2l'^d'd\ (26) 

p'' = -o^' = -o^m' - orU^'-dr, (27) 

^UlU'-'dr, (28) 

which, together with the first of Eqs. |2Jl, yield the following constitutive relations in momentum space 



D = (^ao + a2k^e j E - (^ctq + mit^^j B + [ioxi + a2c^r j (k x B) , (29) 

H = (^/3o - iaxi^t) B - (k x B) + ctqE. (30) 

In the approximation to order ^ here considered, we have ^^fc^ ~ ^iJ^ and thus we can write 

D = rfi(tj)E + id2(t^)B + d3(cj)|(k x B) , (31) 

H = /ii(w)B + i/i2(w)E + i/i3(w)C (k x B) , (32) 

where the functions di{u;) and hi[ijj) depend only on oj and admit a series expansion in powers of S^uj, characterizing 
each specific model. From Eqs. pi2|l we get the equations for E and B 

idi(k-E) = 47rp(cj,k), (33) 

iwdiE + {h^ik^ - g{uj)) SB + (wdaC+^n) («k x B) = 47rj {uj, k) , (34) 

where we denote 

{d2 + h2)uj = g{uj)l (35) 

The expressions (|29I30|I indeed indicate that the above combination is of order ^. Thus we see that in fact there are 
only three independent functions of lo and k which determine the dynamics, which are 

P = rfi , Q = hi+ ujdjl R = {h^k^ ~ g{Lu)) |. (36) 

Using the homogeneous equation ojB = k x E, which yields a; (k x B) = (k • E) k— fc^E, and charge conservation, 
Lup — k • J = 0, we can decouple the equations for the fields E and B obtaining 

k.E=i^(k.j), (37) 

i (cj^F - fc2g) E + i? (k X E) = 47rcj j-^(k.j)k , (38) 

and 

k • B = 0, (39) 

i {uj^P - k^Q) B + i? (k X B) 47rfc k X j. (40) 

Finally, we can introduce the standard potentials <f> and A 

B = ik X A, E = iwA-ik$. (41) 



5 



As in the usual case we can use the radiation gauge, k • A =0, in which case we have 



[k^Q-u^P)K + R (ik X A) = 47r j- j • k k = 47rj 



(42) 
(43) 



from Eqs. H33I34|I . The presence of birefringence depends on the parity violating term proportional to kR, where we 
have a crucial factor of i. This makes it clear that a possible diagonalization can be obtained by using a complex 
basis, which is precisely the circular polarization basis. In fact, decomposing the vector potential and the current in 
such a basis 



and recalling the basic properties k x A^ 



A = A+ + A , Jt = +jT, (44) 
= — iA+, k X A^ = iA^, we can separate (|43ll into the uncoupled equations 



[fc2Q-w2p + Afci?] A^ = 47rj^, A = ±l. (45) 
In terms of the basic functions introduced in the constitutive relations (|31I32|I , the factor in (|45|l is rewritten as 

k'^Q - u^P + XkR^ A/i3fc3|+ (^hi + d3Ujt) - Aff^A: - diuj^. (46) 

This is the key expression to obtain the Green functions and the refraction indices. 



III. THE ENERGY-MOMENTUM DENSITY 



In order to derive the Poynting theorem for this generalized electrodynamics we still need to make explicit the 
remaining dependence upon dt of the operators ao: /^O: ""Oj Q^i, Pi, cri, and a2, in Eqs. (|26I28|) . We denote them 
generically by Co, Ci ^'^d C2; respectively. Following the symmetry conditions imposed by Ea.()13|l we have, up to 
second order in ^, 



Co — Coo + C02C dt^ Ci — Cio, C2 — C20, 
where Coo, C02, Cio and C20 sue now constant coefficients. Here we take 

«oo = PoQ = 1, Coo = 0, 



(47) 



(48) 



to recover the usual vacuum as the background for ^ = 0. A long but straightforward calculation starting from 
Maxwell equations (|1I2|I in coordinate space leads to energy conservation in the form 



ftu + V • S = -j • E. 



(49) 



Outside the sources, where we have the additional property ^^V • E = 0, the corresponding energy density u and 
Poynting vector S are 



, = i (e2 + b2) - |laioE • V X E - f^AoB • V x B 



E.a2E-i(atE) 



(50) 



E X B 



^faioE X 5tE 



HP- 
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-dtB X B - E X (V X B) 



aoaC'E X + p^^^-^E x d'^B. 



(51) 



In our case neither the energy density, nor the Poynting vector retain their standard forms (E • D + B • H)/2 and 
E X H respectively, because we are dealing with a higher order theory. This can be seen from the application of the 
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apropriately extended Noether theorem to the Lagrangian density 112() . which includes up to third derivatives in the 
photon field A^. To this end we introduce 

as a more compact way of writting the parameterization H19I20|I , (|21I24() and (|47|l of the generalized susceptibility 
operator. The coefficients -y^^'^^^---^^"^^ ^ A = 0,1, 2, are constants and the gauge invariant energy momentum tensor is 



ope r 

m 



This satisfies (?rTJ — outside the sources, in virtue of Maxwell equations The standard contributions arise 
precisely from the first two terms in the right hand side of (|53|l . In fact we have 

-L - Fo„i/°" - i (E • D + B • H) , (54) 

-FoaH'"' = E X H. (55) 

It is convenient to remember that these models have to be interpreted as effective theories, in which the dynamics of 
quantum gravity fluctuations are partially incorporated via the modifled constitutive relations (|3I4|I . In this sense they 
are anologous to the Euler-Heisenberg effective Lagrangian for the electromagnetic field including quantum corrections 
through non- linear terms, or to the purely classical Lorentz-Dirac dynamics for a charge including the radiated field 
through a third order derivative term. Thus the density of energy and momentum given by the Noether theorem, u 
and S of Eqs. H50I51|I in coordinate space, contain in principle the contributions of the electromagnetic field plus the 
quantum gravity fiuctuations. These magnitudes satisfy the continuity equation l|49(l . Hence , if we consider a volume 
V outside the sources such that the total flux of S in its boundary is null, it results that 

U = I (fxu (56) 
Jv 

is a constant of motion despite damping factors appearing in a plane wave field. The density u does not correspond 
to the usual expression in electrodynamics, {E'^ + B'^)/2, which does not lead here to a constant of motion. In terms 
of the analogy with an electromagnetic field propagating in a medium, this can be interpreted as the typical effect 
of an active medium that interchanges energy with the field. If this active medium absorbs or cedes energy to the 
electromagnetic field it is a matter of the second law of thermodynamics, which is not explicitely implemented in 
these models. 

We will demand that the physical quantities u and S be causally related to the sources. The fields E and B are 
causal by their construction in terms of retarded Green functions. Since we are dealing with an effective theory with 
an energy cutoff, neither the time nor the spatial derivatives of a causal function are causal, as discussed in Appendix 
A. In this way, to meet the above requirement we have to impose causality, via the Kramers-Kroning relations, to a 
sequence of derivatives of E and B dictated by the order in ^ to which we want the theory to hold. We have explored 
this situation only to first order in ^, where we can show that it is enough to demand also that D and H be causal in 
order to guarantee the causality of the energy density and the Poynting vector. Keeping the approximation to first 
order we can use the zeroth order relations 

V • E = 0, V • B = 0, (57) 
V X E = -dtB; V X B = ^jE, (58) 

in all derivative terms containing a factor ^. In this way we can rewrite u and S only in terms of V x E and V x B. 
On the other hand, from the constitutive relations 

D = E + CTioC'V X B - aio^V x E, (59) 



H = B - /3io^V X B - CTio^V X E (60) 

m express the energy density and Poynti 
and H. The causality requirement for D and H are explored in the next section. 



and assuming aio/3io + (cio)^ Oi can express the energy density and Poynting vector only in terms of E, B, D 
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IV. REALITY AND CAUSALITY CONSTRAINTS 



The reality constraints are straightforward. From the constitutive relations H31I32|I the fields D and H are real for 
E and B real if 

dl{uj) = di(-w), d*2{uj) = -d2{-uj) , dl{uj) = -dsi-Lo), (61) 
hl{uj) = hi{~uj), hl{uj) ^ -h2{-uj), hliu) ^ h3{-uj). (62) 

Even though wc will subsequently extend the variable u to the complex plane, the relations (|61I62I) are written for 
real lo, which represents the physical frequency. Nevertheless, in the case of an absorptive background with a real 
frequency, the wave vector k(ijj), and hence k{oj) = -\/k • k together with the index of refraction n(w) = fc(w)/a;, will 
be complex in general p^ . Unless explicitly stated, all our subsequent expressions involving the frequency are written 
for real lo. 

To have a causal model it is necessary that not only the propagation of the vector potential A, together with E and 
B, is causal, but also the response of the medium coded in the constitutive relations. In the following we discuss this 
last point, finding the causal constraints which must satisfy the functions that appear in the constitutive equations. 
To this end we generically denote the functions di(oj), «d2(w), uj£,d3{oj), /ii(a;), i/i2(w), and iwS,h^{Lo) by ^{u!). In this 
way the reality conditions are 

Rei?(w) = Rez9(-w), Imi?(tj) -lmd{~uo). (63) 

A general discussion of causality and dispersion relations is given in Ref . [T9l| . In brief, the causal response of the 
effective medium depends on the absence of poles of ^{lo) — Rei?([j) + ilmi9(cj) in the upper half complex plane uj. 
In addition, here it is necessary to take into account that we are discussing effective models which lose their physical 
meaning at high frequencies. Therefore, we must consider a finite range of frequencies instead of the usual infinite 
one, and thus causality is characterized by the generalized susceptibility theorem [T ^ [l0 | 

Re.(.)-. = ^p/ d.'r:i:^^^P (64) 









dLu' - 




-n 








/ duj' 


TT 


J-n 



UJ' — U) TT 







T ai \ ^ I 1 ,Rei9(w )-K A 

Iradluj) = 1 duj' h— , (65) 

~ ' UJ'—UJ UJ 

where the cutoff il defines the region where the effective theory holds, < fl, k — limi^^o Re i5(a;) and A = 
—i lm\j^() [uj'd{uj)]. We assume that the contribution of the path closing the loop in the upper half complex plane lo is 
negligible. We have evaluated the above dispersion relations taking an inverse power law behavior of i?(w) for \lo\ > 
and matching the functions at lo = ±f2, obtaining response coefficients which are of the same order of magnitude as 
those in the proposed simpler approximation. Via the dispersion relations, the introduction of a cutoff in lo leads to 
corresponding cutoff in k{Lo). This cutoff induces an effective granularity in space-time at small distances, and can be 
readily seen as follows. Lets us consider a given function in the fc-space, f{k). If there is no cutoff in this space we 
have the corresponding Fourier transform to the coordinate space 

1 r°° 

fix) = — dk e-^'^^fik). (66) 
But if there is a cutoff x in the /c-space , i.e. — x < fc < x, instead of this we have 

fix) = r dk e^'^fik), (67) 
which, using Eq. H66|l . becomes related to fix) by an integral transform 

m = r dx'fix') ''^^^-^'\\ (68) 

J-oo TTix-x') 

In the integrand of the right hand side of Eq. 1681) there now appears a smeared distribution around x, with a width 
X^^, which can be interpreted as an effective granularity in the coordinate space. This space-time granularity is telling 
us that we can not probe space-time at distances smaller than 1/x through the effective theory. 
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If we restrict ourselves to perturbative expressions, we can establish the causal structure of the coefficients in terms 
of an expansion in ^. To be precise, let us write the real and imaginary parts of -d as 



Imi? (w) = 0iuj + 6*3^3 + 0(050;^). 



(69) 
(70) 



In terms of the scale ^ the notation is 9k — i^fcC*^- Moreover, i?(ijj) has no pole at a; = 0, hence A — 0, and we have 
Rei9 (fl) — K a.t ui = U. Thus we get 



02^^ + -e^n* 



n 



This allows us to identify 



ei = - 



and solving for 64 we can predict the third order contribution to 1ijii3{uj) to be 



h + ^02n 



Thus the final expression up to third order in ^ is 



^1 , 161^2 ^,.P3 

UJ I; . 



This is the structure imposed by causality on every coefficient in the relation of E, B with D, H. 



(71) 



(72) 



(73) 



(74) 



GREEN FUNCTIONS IN THE RADIATION GAUGE 



The retarded Green function for the potential A and the field strengths E, B (sec Eqs. (|38I4UI43|) ') can be written 
in terms of the circular polarization basis 



GlfM = J 0e-Grf (c.,k) = i J ^e-- ^G^.,k) (^..-^ + *A.,.,| ) , (75) 
where ki = fci/|k|, k = |k|, and G^{LLj,'k) is obtained from Ea. l45|l 



To obtain the causal Green functions the analytic continuation uj uj -\- ie is taken. Only the poles in the upper half 
plane of the variable k make a contribution to the integration. The first step is to rewrite the denominator in Eq. 
()76|l in a convenient form. This is done by successive rescalings leading to 

Qk^ - Puj^ + XkR = --nluj^Qa{M^ - Ma){M^ - M+)(M^ - M_), (77) 
3 

where we introduce the notation 

Q = hi+d3uj^, X^^ujQ"^, nl=diQ^^, 

a = ihsnox, = Xk (riQu;)^^ . (78) 

The corresponding roots are given by 

1 1 1 + ac A 



24/3a V J 22/3a^' 



+ 
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with 



2 + 3ac - 3a^ + J (2 + 3ac - ?,a?Y - 4(1 + acf 



1/3 



(80) 



and c = gx (^^^o) To study the modifications to the dynamics it is enough to expand each root in powers of the 
small parameter x 



Mo ~ jx'\ M±~± 1 + i (/3i - ai) (^Xx + j (5/3i - &i 



X 



(81) 



where /3i = h^UQ and ai = g/(uj^no). Since the parameter A and the momentum k appear only in the combination 
\k it is clear that we have the symmetry property 



G'-'(w,fc) = G--'(cj,-fc). 

This relation will be useful in the final calculation of the Green functions G'^{lo,y). 
The integral in H75|l can be written 

Glf (^, r) = [G^f ] ^ (c., r) + [G^] ^ (c., r) + [G^] ^ {cu, r). 

Each vector kj can be obtained by inserting —idj outside the integral. In this way we have 

e^^-^Y.G\u,k) = ——^5,u / kdk{e^^^) [G+(c., fc) + G-(c., fc)] 

/TT , ' 27r) J-oo 



(82) 



(83) 



2[G:t%{u:,v)^d,dk 
2[G:^\{uj,v)^e„A 



(Pk e''^" 
(Pk e'^" 



i 1 



■i-n^nk) / kdke'''''\G+(uj,k) + G-(Lu,k)] 

r {2nY 

^^G^iuj.k) ^ --^-^it^jkUj I kdk{e'^'') {G+iuj,k) ~ G-{LU,k)) 



i2nr 



where n.^ — Xi/r. Thus we arrive at 

(^,r)=- 



i 1 



(27r) 



1 f°° 

o i^tk - niUk + iXejpknp) / kdke'''^'' G^{uj,k), 



which allows us to identify G^{u},v) 



G^{u,v) = ---^ 
r{2^f 



kdke'"'' G''{uj,k). 



(84) 



(85) 



(86) 



Let us emphasize that the factor e''^'' forces us to close the integration contour by a circle at infinite in the upper half 
complex plane, thus picking up the poles in this region. Our next step is to perform the integrals in To this end 
we recall that 



G^{Lj,k) = -- 



1 



iluj^Qa {M^ - Mo) (M^ - M+) [M^ - M_) ' 



(87) 



Our description is an effective one valid only for momenta k << Mp. According to Eqs. (|81|l . the pole at AIq 
corresponds to the momentum value 

\ko\^\Qh^'\r'- (88) 

In the approximation considered here this pole can be taken at infinity and its contribution to the integral can be 
neglected. The two remaining poles, which are the ones that contribute to the integral, are located at very small 
displacements with respect to \k±\ ~ uquj << Mp. In this way we get 



G^{uj,k) = 



1 



iluj'^QaMo (M^ - M+) {M^ - M_) ' 



(89) 
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Let us consider first the case of G"*", where we have 

G+{uj,k)= ^ 



1 



QuMq (fc — a;rioAf_|_) (fc — ujnoAI^) 



(90) 



From the leading order expressions in (|81|l we conclude that the pole that contributes in this case is {uj + ie) noAf+. 
Thus the resulting integral is 



G+(tJ,r) 



1 



inr QaMo {noM+ ~ n^MJ) 



(91) 



The second integral proceeds in an analogous way, but now {— (uj ^ it) n^M J) is the pole that contributes. It is 
convenient to define the refraction indices in the circular polarization basis as 



(92) 



The minus sign in n_ is due to the fact that starts with a —1. Up to the order considered, they are given by the 
explicit expressions 



nx{Lu) = no 



1 + A /3i - ai 



I + (/3i - 5i) (5/3: 



1 - ai 1 ^ 



From here we can obtain an explicit second order power expansion in ^ 



where 



(93) 



(94) 



'^1 = 77 ("10 - Pw) 



"-2=0 [("10 - Pw) ("10 - 5/3io) + 4 (/3o2 - "02) 



(95) 
(96) 



With this notation we finally get 



2nx 



AtttQ (n_ + n+) 



(97) 



where we have considered that the dominant term in Mq yields ail/o = 1. The reality conditions Hfill62|l impose the 
relations 



[G^(c^,r)] =G-^(-c.,r), [n^iut = n^x{-uj). 



(98) 



The low energy effective character of the theory implies the existence of a frequency cutoff Q,. Taking this into account, 
the space-time Green function is 



1 



n 



2nx 



1 

47rr 



^ Q{n-+n+) 
n 

'1 



duj 



n 



+ AniLj|- {a2o - P02) (^^tj 



iuj 1 1+Ani LJ^+n2 (cj^^ r 



(99) 



where t = t — t' . The choice of the poles in the complex plane uj is consistent with the causal behavior of the Green 
function. But the frequency cutoff could introduce some violation of causality. To investigate this possibility, we can 
compute the time dependent Green function by expanding the integrand in powers of ^ 



G^(r,r) 



1 

47rr 
1 

27rr 



du! <1 + Xfii (1 + iusr) Lu£^ 



"20 - Po2-i ("2 + nf) ruj + ^n^r^w^ 



1 - iXnii,{l + rdr)dr + l'^ 



"20 - P02 - {n2 + n\) rdr - -;^n\r'^dl 



dt 



2I sin (r — r) r2 



(100) 
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This expression shows that the main effect of the cutoff is to spread the propagating field around the hght cone, 
within a wedge defined by r ~ r ± ti/2Q,. Returning to G^{u>, r), Eq. (|97|) . we can characterize the effect of the cutoff 
in the causal behavior of the Green function using the generalized susceptibility theorem flOj, a generalization of the 
Kramers-Kronig relations. In terms of the frequency lo, the real and imaginary parts of the Green function are 



ReG\LO,r)^ — 



Im G^{Lj,r) ~ 



1 + ( ^20 - Po2 - ("2 + nl) rdr - -nlr^dl ] 8^ 



coswr + Ani^ (1 + rdr) dr sin (wr) j- , 

(101) 



1 + ^20 - /3o2 - («2 + nj) rdr - ^nydf. dl 



sinwr — Ani^ (1 + rdr) dr cosLor 



(102) 



According to the generalized susceptibility theorem the Green function will be causal if 

,ReG(tj')-ReG(f7) 



1 r 

Im G{uo) = P / duo 



(103) 



or more explicitly when 
ImG^(cj) = 



1 + "20 - /3o2 - (712 + nl) rdr - T^n^dl d, 



P 



duj 



I cos {iij'r) — cos (i7r) 



^ ^ \ ^ f^^ , , sin (cj'r) — sin (fir) 
+Ani.^ 1 + rdr) 9,. P / duj' ^ ^ ^ — ^ 

' OJ—LU 



(104) 



For oj/il <^ 1 the integrals reduce to 



P 



, cos (uj'r) — cos (fir) 



duj 



2 (1 — coscji?) — cos (fir) — 



TT + [{fir) (cos ilr) 



smujr, 



„ , , , sin (w'r) — sin (fir) , ^ uj ^ 

P I duo' ~ 2 (1 - sinwr) - cos (fir) + 

Q uo' — UJ fi 



TT + [(fir) (cos fir) — sin fir] y ^ 



(105) 



Furthermore, in the case of a radiation field fir ^ 1 and hence the factors cos fir and sin fir become strongly oscillating, 
nullifying the contributions of the terms where they appear (which also have a factor (cj/fi)", with n>\). Thus we 
can take 



P I duo 

-a 



,cos (oj'r) ~ cos (fir) 



„ , , .sin (oj'r) — sin (fir) 

P duo' i '- i — ^-TTCOswr. (106) 

l-n Lo' -UJ 



Replacing these integrals in (|l()4|l . we finally get that if Re G^{uj) is given by Eq. (|l()l|l . the imaginary part of the 
Green function, Im G^{uo), must be 



Im G^(cj) 



47rr 



1 + "20 - /3o2 - (n2 + nl) rdr - ^nyd^ d, 



sincjr — Ani^ (1 + rdr) dr cosuor > , (107) 



which coincides with Eq. (|102|l , obtained by direct computation. This result shows that the cutoff does not introduce 
any significant causality violation provided that fir >> 1 and uo/fl << 1. 



VI. FINAL COMMENTS 



In the preceding sections we have discussed the implications of reality and causality of the electromagnetic fields 
upon the structure of the constitutive relations and the Green functions in quantum gravity inspired modified elec- 
trodynamics. To illustrate these implications we now apply our results to some particular models already found in 
the literature. 
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Let us start with the Gambini-PulHn electrodynamics [20j, where the constitutive relations are 

D = E- 2i^wB + 4^2^ k X B, H = B + 2i^kxB. 
Introducing the dominant contributions to satisfy causality they become 



D = E - 2i^ujB + 4(1 
which leads to the refraction index 



16 iuj 



C^cj k X B, H = B + 2i^k X B, 



o ~o 64 w „ ~o 
2Acj^ + 4tj2^2 _ j_ 2^2 

3tt il 



(108) 



(109) 



(110) 



indicating that absorption is generated at the order {^lu^'^ 
The second case is the model of Myers and Pospelo 

D = E + i^LuB 

Here the constitutive relations become 



which yields 
H = B + i^ujE. 



D = E 



1 



H = B + 2 1 + U^^E 



after the causality requirement is imposed, and the resultant index of refraction is 



n\{LL>) 



1- A 1 



UJ 



UJ 



2c2 



(111) 



(112) 



(113) 



In this case the absorptive terms begin at least to order uj^ and causality only introduces fJ-dependent corrections to 
the velocity of propagation. The two models discussed previously present birefringence. 

The third model that we consider is the flat space version of that of Ellis et al. 0,13,13, which has been used to study 
in a rough way the phenomenology of the full underlying model, leading to robust limits on Lorentz violation. From 
the theoretical model of a recoiling D-particle in a quantum gravitational foam, a simple effective model formally 
analogous to electrodynamics in a medium has been heuristically introduced in Ref. . Such an approximation has 
been used in conjunction with observations of gamma ray bursts to set limits on the quantum gravity scale Afj^. 
This effective model is characterized by constitutive relations Q 



D = ^ + H X G, 



B = ^ E X G, 



(114) 



where 1/Vh plays the role of the electric and magnetic permeability and is taken equal to one to have the same 
permeability as the classical vacuum. In references ^ 0| the vector G, which originally represents the recoil velocity 
of the I?-particle (being proportional to the momentum transfer over M) is very crudely identified with G ^ k/i\f , 
where k is the momentum of the propagating photon and M is the effective mass of the D-particlc in the quantum 
space-time foam. In our notation this amounts to require 



G=/k 



and the above constitutive relations turn out to be 



D = E+(/2.^-/)(kxB), H=(l-/c^)B. 
To lowest order, the reality conditions for the fields imply that we must choose 

/ = 

with a being a constant real number. Imposing causality we get 



D = E 



^(kx B) 



H 



B. 



(115) 
(116) 

(117) 
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The dominant contributions to the refraction index do not depend on the causal corrections and are 

nx{Lo) = l + iaLui-a^LO^(_'^. (118) 

Contrary to the previous cases, the above indicates that absorption is already present to linear order (a point previously 
missed) and that the propagation properties of this effective model exhibit corrections starting from second order in 
^. In fact the photon group velocity is given by 

c(^) ^ ^IIM = 1 _ 3{a^Lu f, (119) 

and not by c{uj) — co(l — uj/M), with M real, as was originally assumed in setting the observational bounds 

To summarize, our discussion is based on the consideration that these modified electrodynamics are effective theories 
with physical meaning, and hence that in particular the electromagnetic field strengths E and B, and the excitations 
D and H, must be real and its propagation causal. We have shown that these requirements impose constraints on 
the structure of the constitutive relations and lead to additional contributions to the photon propagation velocity 
and, much more interesting, to absorptive effects in the propagation of the fields. The contributions originating from 
causality depend on the range of validity of the effective theory, coded here by the cutoff f2. 
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Appendix A 

In this Appendix we analyze the effect of a frequency cutoff on the causal character of the derivative of a causal 
function f{t), defined by a generalized susceptibility {}{t) 



f{t)^ I dTd{t-T)g{T). (120) 

The derivative is 

dtf(t) = / dr dtd{t - T)g (r) , (121) 



and therefore the problem reduces to studying the causal character of the derivative of '!?(i), considered itself as a 
susceptibility function. We are dealing with real functions, so we will consider here a real susceptibility , d{t) = ^}*{t). 
Introducing the Fourier transform 



= / dt e*"*i9(t), (122) 

— oo 

we get in this case 

Rei^(w) = Rei9(-w), Imi?(tj) = -lm§{-Lj). (123) 

If ^{t) is causal, its Fourier transform satisfies the Kramers-Kronig relations, which in the presence of a frequency 
cutoff become 

1 T™ „<n, ,'\ 

K, (124) 





•n 




duj' - 




-n 








/ du' 


n 


J-n 



— UJ 



. , , Rci?(cj') - K A , ^ 

Imi9w= P\ dJ + . 125 
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where k = Rez?(n) and it is assumed that = i-^ at ~ 0. The Fourier transform of the derivative of the 

susceptibility function is given by 

d{uj) = (126) 

and hence 

Re{){uj) =ujlmi9{uj), (127) 
lm{){uj) = -ujRei}{uj). (128) 

In order for 'i?(tt') to also be causal, it must satisfy the corresponding Kramers-Kronig relations 

k, (129) 

(130) 









duj' - 




-n 






-Ip 


/ dijj' 


TT 


J-n 



LU' — LU 

, Rez?(cj') -k A 



where k — Rei9(r2) = u!lm'd{Lu)\^^ and we are assuming that has a pole of order n at uj = 0. Rewriting Rei9(cL;) 
and lm{){u!) in terms of Re-d{u!) and lini3{uj), according to (|127|l and (|128|l . these last equations yield 

Imi9(w) = P\ duj' — — + -, (131) 



TTUJ J _Q UJ' — UJ UJ 



Re^(c.) = —P / du^' —. (132) 

TTU) J _Q Uj' ~ UJ Ci;"+-^ 

Finally, comparing Eqs. p24ll25() and p31ll32() . and considering that Im 19(0-1) is an odd function, we get 



A K , /I - uj/Q 
In — - 

TTUJ \1 + UJ/Q 



1 , , ,^ K , fi-uj/n\ A k , , 

— duj'Rediuj') = —lni ^ - + -. (134) 

TrUJJ_Q TT \l+Uj/ltJ UJ"^ UJ 

The first relation implies 

A = K = k^O, (135) 

and thus the second one leads to 

/n 
duj' Re^iuj') = -ttA. (136) 
-n 

When there is no cutoff at a finite frequency the factor In ^ becomes 0, and the constraints are 

K = i = 0, (137) 



d^' Re ,?(.;') = -(--^) 



(138) 



For example, let us consider what happens when i!){t) is the causal Heaviside function, such that ^{uj) ^ i/uj and 
liiuj) = 1. In both cases, H. ^ 00 or fl finite, we get 

K = A = Q A = k=l, / duj' Re d{uj') = 0. (139) 

J-n 

When — 00 the constraints (|137I138|I hold and are indeed satisfied, and hence dt9{t) — S{t) is also a causal 
susceptibility. But if there is a finite frequency cutoff the corresponding constraints are given by p35ll36|) , which 
are not satisfied, and S{t) becomes a non causal susceptibility. 
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In the case of a field theory we have spatial derivatives together with the time derivative. But the former derivatives 
are related to the latter by the equations of motion that drive the evolution of the system. Thus, to discuss causal 
characteristics it is necessary to take into account both types of derivatives, linked at the Fourier transform level by 
the refraction index that relates frequency and momentum. 
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